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ABSTRACT 



Aircraft sctructures may be modelled by lumping the masses at particular strategic 
points and the flexibility or stiffness of the structure is obtained with reference to these 
points. Equivalent moments of inertia for the section at these positions are determined. The 
lumped masses are calculated based on the assumption that each point will represent the 
mass spread on one half of the space on each side. Then these parameters are used in the 
diferential equation of motion and the eigen characteristics are determined. A comparison 
will be made with results obtained by other established methods. 

The lumped mass approach in the dynamic analysis of complicated structures 
provides an easier means of predicting the dynamic characteristics of these structures. It 
involves less computor time and avoids computational errors that are inherent to the 
numerical solution of complicated systems. 

INTRODUCTION 

The mass of the hypersonic plane is continuously distributed over the entire 
structure. Consequently, the real structure has an infinite number of degrees of freedom as 
far as the dynamic behavior is concerned. However, in the dynamic analysis of structures, 
it is possible to replace the real structure with an ideal one consisting of a number of 
lumped masses. These are assumed to be connected to one another through elastic massless 
elements which, to a certain extant , retain the actual behavior of the original structure. The 
method of idealizing actual structures bears significantly on the final results in any vibration 
analysis, and the selection of the method and the number of lumped masses for the system 
has to made while taking into consideration the various aspects of the structure under 
study. The skill and experience of the analyst are very helpful in obtaining the best ideal 
model for the structure. 

In idealizing the hypersonic plane, there are certain assumptions which have to be 

made: 

a. Idealization of real strucutres is limited to those structures which deflect in a 
linearly elastic manner. It is possible to extend the procedure to structures 
loaded in the plastic region, but the solution of such structures is more 
complex. 

b. To be idealized, a structure must be stable under static loads. This condition 
applies for both determinate and indeterminate structures. 


m 


c. All structures demonstrate a certain amount of damping when they are 
subjected to dynamic loading conditions. Such damping in structures is 
controlled by structural hysteresis and by external friction. In the dynamic 
analysis of ordinary structures damping may be nec glee ted in determining 
the natural frequencies, but it must be included in the evaluation of mode 
shapes under resonant conditions. 

DEVELOPMENT OF THE METHOD 

Consider that the frame of the hypersonic plane is fixed at the narrow end to act as a 
cantilever and that the masses are lumped as seen in Figure 2. In this study the shear and 
rotary inertia effects are ignored, the dynamic loading on the cantilever beam is the inertia 
of the moving bodies. The inertia force due to each body is expressed as 

Mi Zim - - Mj CO^Zim ( 1 ) 
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Figure 2 

where: 

Mi = the lumped mass at a point i 
Zjm = the deflection of point i in the mth mode 

Zim = the acceleration of point i in the mth mode 
0 = the circular frequency of the system vibrating in the mth mode 

The deflection equation for the structure under dynamic loading due to inertial 
forces is expressed in the following form: 

{Z m }=co2 [A][M]{Zm} (: 


{Zm} = a column matrix of the displacement of the structure in the mth mode. 

[A ] = a square matrix of the flexibility coefficients of the structure 
[ M ] = a diagonal matrix of the mass of the structure. 

Equation ( 2 ) can be expressed in the alternate form as in 

0 = [ D ] {Z m } (3) 
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where 


[DJ =[[I]-g£[A][M] (4) 

To obtain a non-trivial solution for Equation ( 3 ), the determinant of matrix [ D ] 
must be identical to zero. 

0 = I D I (5) 

The expansion of Equation ( 4 ) yields the characteristic equation for the stucture 
which is a polynomial. The nth degree of the characteristic equation is equal to the rank of 
the matrix [ D ]. The roots of this equation represent the eigen values of the structure. 

DEVELOPMENT OF THE FLEXIBILITY MATRIX [A] 

The model as shown in Figure 2 has displacement in one plane only. The 
displacements are considered to be lateral and rotational ones. The shear displacements are 
neglected in this model. The orientation of the model displacements is shown in Figure 3. 



Figure 3 


SYMBOLES 

N = Number of lumped masses 

L = Length of beam 

I = Moment of inertia 

E = Y oung's modulus of elasticity 

mj = ith mass of the beam/mass moment of inertia 
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p - Density of the beam per unit length 

Considering that the beam is divided into equal segments along the longitudinal 
axis, the length of each segment is given by 

( 6 ) 


The mass of a segment of the beam at any point i is 

2.788 r , , , 24 

+ m ( H4 - Xi)] , i = 1, N 

The function for Xj is given by the following: 

Xl = 2N 

L . „ xr 

Xi - Xj_i + ^ , i — 2, N 


(7a) 

(7b) 

(7c) 


The mass moment of inertia is assumed to be that of a bar with a uniform mass over the 

length of the segment. It is given by the expression 

l 2 

mmii = mi] j2 . i = N + 1 , 2N ( 8 ) 

where 1 = ^ , mij = raji = 0 and mmij = mmji = 0 


The flexibility matrix [A] in Equation 2 is expressed as follows: 


[A1 = Aj| Ak ( 9 ) 

where 

[ An ] = is NxN matrix that represents the translational diasplacements 
due to unit lateral forces. 

[ An ] = is NxN matrix that represents the rotational displacements due 
to unit lateral forces. 

[ A 21 ] = is NxN matrix that represents the translational displacements 
due to unit rotanational forces. 

[ A 22 ] = is NxN matrix that represents the rotational displacements due 
to unit rotational forces. 

The matrix [ An ] is generated from the following equation: 

a U “ M 1 4 Fp [ ( 2N + 1 - 2i > 3 - 3(2N + 1-21)2 (j-i) + 4 0 -i) 3 ] ( 10 ) 

a Ji = a U (11) 

for i = 1 , N and j = i, N. 

Matrix [ A 12 ] is obtained by taking the derivative of Equation (9). 

k>j = • 2Ei ( 4 nI (2N + 1 - 2i) 2 - 4(j -i) 2 ] ( 12 ) 
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for i = 1, N and j = i , N. 

bji = _b« , i = 1, N and j = 1, i ( 13 ) 

The development of the matrix [ A21 ] follows from the application of unit moments at the position 
of the lumped masses on the beam and finding the lateral displacements that ensue from these 
actions. These displacements are given by 

C « = 2EI ^ 4N2 )[ ( 2N + 1 - 2i ) - 2(j - i)] 2 ( 14 ) 

for i = 1, N and j = i, N. 

. . L j 

(15) 


Cji - Cii + (i - j ) ^ d ;i 


for i = 1, Nand j = i, i - 1 

The roatational displacements matrix [ A22 ] due to rotational forces is expressed as 

dy = gj ( 2N ) t ( 2N + 1 - 2i ) - 2( j-i )] (16) 

for i = 1, Nand j = i, N. 

dji = dy , i = 1, N and j = 1, i. ( 17 ) 

Equations ( 10 ) through ( 17 ) define the flexibilty matrix for the entire structure. Having 
determined matrices [ A ] and [ M ], substitute them in Equation 4 and solve for the eigen values. 

EXAMPLE 


Consider that the hypersonic plane model is subdivided into three parts. The [ A ] and [ M ] 
matrices are given below. 

An A12 

A 2 i A 2 2 
125 54 7 

n a a 



[A] = 


W&==z 

m 

[ An ] = 

L 3 

648EI 

m 

i 

[A 12 ] = 

L 2 
‘ 72EI 

E3 

[ a 21 ] = 

L 2 

72EI 


[ A 22 ] = 

L 

6 EI 


25 21 9 
9 9 5 


25 9 1 
21 9 1 
9 5 1 

5 3 1 
3 3 1 
1 1 1 
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where 


[M] = 


Ml 0 
0 M2 





K4 

Ej 

m 


I 

I M 


and 


[M t ] = 


1 

386.4 


33.0 0.0 
0.0 25.3 
0.0 0.0 


0.0 

0.0 

17.7 


[M 2 ] = 


1083 

386.4 


33.0 0.0 
0.0 25.3 
0.0 0.0 


0.0 

0.0 

17.7 


The modulus of elasticity E = 10 7 psi. The equivalent average moment of inertia I = 42.5 
in 4 . The length of the beam L = 1 14 in. With the values of the matrices [ A ] and [ M ] known, the 
solution of Equation ( 5 ) provides the following frequency results. 

The natural frequencies ensuing from the model with three lumped masses are given in the 
first row for the first six modes. Values in subsequent rows correspond to models of 4, 5, 6, and 7 
lumped masses. 


Natural Frequency in Hz for a Lumped Mass system 


Number 

Number of Modes 

III 

Elements 

1 

2 

3 

4 

5 

6 

3 

19.54 

146.08 

396.08 

560.60 

1329.9 

2125.6 

4 

19.21 

138.15 

41 1.18 

454.22 

1041.8 

1492.6 

5 

19.08 

134.70 

399.10 

509.07 

779.19 

1661.1 

6 

19.00 

132.87 

391.36 

558.29 

778.32 

1247.7 

7 

18.95 

131.78 

386.46 

603.53 

770.65 

1273.1 


Natural Frequency in Hz for a Finite Element Model 


Number 

of 

Elements 

Number of Modes 

1 

2 

3 

4 

5 

6 

3 

18.11 

132.15 

385.01 

874.00 

1656.55 

3108.26 
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CONCLUSION 


A lumped mass model for the hypersonic transport airplane has been established. 
Algorithms for the determination of the flexibility matrix [A] and the mass and mass moment of 
inertia matrix [M] have been found. The natural frequencies for a 3-lumped mass system have 
been determined using the lumped mass method and the finite element method. The results from 
the two methods converge in the lower three modes and diverge in the upper three ones. The 
lumped mass system requires less computer time than the finite element model. For models with 
large number of elements, the lumped mass system is more efficient Results from both models 
need to be verified experimentally. 


